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Abstract 

We calculate the friction force between two semi-infinite solids in 
relative parallel motion (velocity V), and separated by a vacuum gap 
of width d. The friction force result from coupling via a fluctuating 
electromagnetic field, and can be considered as the dissipative part 
of the van der Waals interaction. We consider the dependence of the 
friction force on the temperature T, and present a detailed discussion 
of the limiting cases of small and large V and d. 

1 Introduction 

Because of its great practical importance and because of the development of 
new experimental techniques, sliding friction has become a topic of increasing 
attention pQ. In this paper we consider the friction force between two solids 
in relative motion, separated by a vacuum gap of width d. This "vacuum" 
friction is in most cases of no direct practical importance since the main 
contribution to the friction force when a body is slid on another body comes 
from the area of real (atomic) contact pQ. Thus, the frictional stress between 
two semi-infinite metallic (e.g., copper) bodies, moving parallel to each other 
with the relative velocity V = 1 m/s, and separated by the distance d — 10 A, 
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is only (see Ref. [2] and below) ~ 10~ 6 N/m 2 . This stress is extremely small 
compared with the typical frictional stress ~ 10 8 N/m 2 occurring in the area 
of atomic contact even for (boundary) lubricated surfaces. Nevertheless, vac- 
uum friction is important in some special cases (see Ref. [2]), and determines 
the ultimate limit to which friction can be reduced. Quantum and thermal 
fluctuations of the polarization and the magnetization of solids give rise a 
fluctuating electromagnetic field. For two stationary solids the interaction 
mediated by this field result in the well-known attractive van der Waals force. 
For two solids in relative motion this interaction will also give rise to a friction 
force between the bodies. The static aspect of the van der Waals interac- 
tion is well understood but there are still controversial results concerning the 
dynamical part. Different author's have recently studied the van der Waals 
friction using different approaches, and obtained results which are in sharp 
contradiction to each other. The first calculation of van der Waals friction 
was done by Teodorovitch e 3]. Schaich and Harris [3], and Pendry[5] argue 
that Teodorovitch calculation is in error. For two metallic bodies Schaich 
and Harris found that the friction force is independent of any metal prop- 
erty, in contrast to the results of other authors. The friction forces calculated 
by Levitovjn], PolevoijT] and Mkrtchiau[8 vanish in the non-retarded limit 
(formally obtained when the light velocity c — > oo). This result is very sur- 
prising (and in our opinion incorrect), since neglecting retardation is a good 
approximation at short separations d between the surfaces, in which case one 
would expect the friction force to be particularly large. Even at large sep- 
aration, where the non-retarded contribution is negligible, our result differs 
form those of Levitov, Polevoi and Mkrtchiau. PendryjS] considered only 
the case of zero temperature, and Persson and Zhang[2] the case of small 
sliding velocities, and both groups neglected retardation effects. To clarify 
the situation we present a straightforward calculation of the van der Waals 
friction based on the general theory of the fluctuating field developed by 
Rytovjlj and applied by LifshitzJTU] for studying the conservative part of the 
van der Waals interaction. In this approach the interaction between the two 
bodies is mediated by the fluctuating electromagnetic field which is always 
present in the vicinity of any collection of atoms. Beyond the boundaries of 
a solid this field consist partly of traveling waves and partly of evanescent 
waves which are damped exponentially with the distance away from the sur- 
face of the body. The method we use for calculating the interaction forces 
is quite general, and is applicable to any body at arbitrary temperature. 
It also takes into account retardation effects, which become important for 
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large enough separation between the bodies. A similar approach was used 
by Polevoi [7] but he obtained a nonzero friction force only in the relativistic 
limit, in contrast to the present calculations and the earlier calculations of 
Persson and Zhang[2j, and Pendry[5 . Polevoi did not give enough details of 
his calculation to compare his theory with the present calculation, but we 
believe that he overlooked effects related to the change in the reflectivity of 
electromagnetic waves from moving bodies, which occur even in the nonrela- 
tivistic limit. In the nonretarded limit and for zero temperature the present 
calculation agree with the results of Pendry. Similarly, in the nonretarded 
limit and for low sliding velocities, we agree with the study of Persson and 
Zhang. 

2 Calculation of the fluctuating electromag- 
netic field 

We consider two semi-infinite solids having flat parallel surfaces separated 
by a distance d and moving with velocity V relative to each other, see Fig. 
1. We introduce the two coordinate systems K and K with coordinate axes 
xyz and xyz .In the K system body 1 is at rest while body 2 is moving 
with velocity V along the x— axis ( xy— and x y — planes are in the surface 
of body 1, x— and X £1X6 S have the same direction and the z— and z — axes 
pointing toward body 2). In the K system body 2 is at rest while body 
1 is moving with velocity — V along the x— axis. Following to Lifshitz to 
calculate the fluctuating field in the interior of the two bodies, we shall use 
the general theory which is due to Rytov and is described in detail in his book 
jH]. This method is based on the introduction into the Maxwell equations of 
a "random" field (just as, for example, one introduce a "random" force in the 
theory of Brownian motion). In the K system in a dielectric, nonmagnetic 
body 1 for a monochromatic field (time factor e~ luJt ) these equations are 

V x Ex = i-Bi 

c 

VxBx = -i-e 1 {u)E 1 -i-F 1 , (1) 

c c 

where E and B are the electric and magnetic fields, £\{oj) is the complex 
dielectric constant for body 1, and F is the random field. Accordingly to 
Rytov the correlation function of the latter, determining the average value 
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of the product of components of F at two different points in space, is given 
by the formula 



i i i , 



Fi(x,y,z)F k (x ,y,z 
Ah Q + n(w)) e"(w)M(z - x)S{y - y')S(z - z), (2) 



n(u) 



e hu/k B T _ I ' 



where T is the temperature and e" is the imaginary part of e = e + e" . The 
function F(x, y, z) can be represented in the form of a Fourier integral, which 
can be written for the half space z < in the form 

/+oo 
g(k)e 4q ' r cos k z zd 3 k, (3) 
-oo 

where a two dimensional vector q is parallel to the surface, k 2 = k z + q 2 , and 
r is the radius vector in the x — y plane. For the Fourier components g(k) , 
the correlation function corresponding to the spatial correlation (3) is 

v h(l + n(u))e"(u) 
9l (k,u)gl(k ,co)) = -A2 \'J J Sik S(k - k ) (4) 

In body 1 (z < 0) the fields E and B can be written in the form [TU] 

n+OO 



/ + CO 
{ai (k) cosk z z + zbi (k) sin k z z) e 4q ' r d 3 k 
-oo 

/+oo 
midy^-^d'q, (5) 
-oo 

Bi = - / {([qxaJ+^[e z xbi])cosfe z 2 

J-oo 

+i([q x bi] + A;^[e^ x a x ]) sin /c 2 2;} e l<l ' r d 3 k 

+- / {[q x Ul ] - Sl [e 2 x Ul ]}e^™ 2 rf 2 g, (6) 

Co" J-oo 

where e z is a unit vector in the direction of the z axis, and 



where the sign of the root is to be chosen so that the imaginary part of s will 
be positive. 

The first ter)s in the expressions (5) and (6) represent a solution of the 
inhomogeneous equations (1). Substituting them in the second equation of 
(1) and writing F in the form (3) , one can find the following relations, 
expressing a x and b x in terms of the Fourier components gi of the random 
field 

^■eigi - q(q • gj - k'ig lz e z 



ai 



Ei(k 2 — U0 2 Ei/c 2 ) 



(8) 



k 

hl = -TTT2 V / 2\ Mq ' Si) + Qgu] ■ ( 9 ) 

The second integrals in (5,6) represent the solution of the homogeneous 
equations (1) (i.e. the equations with F omitted), and describe the plane 
wave field reflected from the boundary of body. The condition for transver- 
sality of these waves is 

ui • q-siui z = 0. (10) 

In the space between bodies(vacuum) e = 0, F = and the field in the K 
system is given by the general solution of the homogeneous equations, which 
can be written in the form 



Eq 



/+oo , . 
U(q,to)e tpz + w(q,cu)e- tpz \ e l * r d 2 q, (11) 
-oo ' 



B 3 = - / ([qxv]+p[e |X vr 

+ ([q x w] - p[e z x w]) e^ 2 } e^ r d 2 q, (12) 

where 

P = (13) 
and v and w satisfy the transversality conditions 

v ■ q+pv z = 0, w • q— pw z = 0. (14) 

The boundary conditions on the surfaces of the media are the requirement 
of continuity of the tangential components of E and B in the rest frame of 
respective body. In the K system on the plane z = for a given value of q 



5 



it is convenient to write the corresponding equations for components of the 
fields along the vectors e q = q/q and e n = [e 2 x e q ], this gives the following 
equations 

r+oo 

/ aiqdk z + Mi q = Vq + Wq, 




a ln dk z + u ln = v n + w n , 



+00 

qa lz - k z biq)dk z + qu lz + si^i q 



/+00 
-k z b ln dk z + si^in = -p(u n - w n ), (15) 
-00 

where ai q = e q ■ ai, a in = e n • a! and so on. In what follow we shall need 
only the field between two media. Using the transversality conditions (10) 
and (14) and the expressions (8,9) from the equations (15) we can obtain the 
following equations 



r°° 1 

PSi / (50i*(q,fcz, ^) - sigiq(q,k z , u)) — j' 

J — 00 fc 2 — Si 

- (si + pei) u q (q,w) - (pei - Si) «; q (q,w), (16) 



/w\ 2 p g ln (q,k z ,u) 

= ip + si)v n (q,o;) + (si -p)w„(q,o;). (17) 

In the if' system the Maxwell equations have the same form (1) and in 
the second medium ( the half space z > d), the field E 2 , B 2 is given by the 
same formulas (5-9) with x— coordinate changed to x , the index 1 changed 
to 2, coskzZ, sm.k z z replaced by cos/c 2 (z — d), sinA; z (,2 — d) and change in 
the sign of s ( the "reflected" waves now propagate along the positive z 
direction). In the space between media in the K system the field is given 
by the same formulas (11-13) with x changed to x , v, w replaced by v , w . 
The relations between the field in the K and K' systems are determined by 
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Lorentz transformation. Neglecting by the terms of the order (V/c) 2 these 
relations are given by 



v'(q, J) = v(q, u) + — [e x [k x v(q, u)]] , 

UJ 



V 



(18) 



w'(q',u/)= w(q, u) + — e x kxw(q, a;) 



V 



UJ 



(19) 



where k = (q,p), k = (q, - p), J = ui - q x V, q' = q- (Voj/c 2 ) e x . 

In the K' system the boundary conditions at the surface of body 2 at 
z — d give the equations 



2 dk% 



= (pej - s 2 > ,(qV)e^ + (s 2 +pej)«, ,(q>V ipd , (20) 



2 



roc 



5 2 



^2n'(q ) ji / - \ ' / ' '\ ipd 



+ ( S2 -+p)w' n ,(q,uj')e-^ d , (21) 

where e 2 = 62(00 — q x V), 

s 2 = V^'/c) W) - q' 2 = ^^f^Mw - q*V) - l)+p\ (22) 

p is invariant under Lorentz transformation. Now from the equations (18-19) 
with accuracy to the terms of the first order in V/ c we have 

vjiq ,u ) = (v -e q >) « u q (q,w) + -^-v n (q,u>), (23) 

' / ' \ w uiq v V 
V = (v • e n » -u n - (24) 

The similar equations can be written for w' ,,w' n , .After substituting (23) 
and (24) into the equations (16-17) and (20-21) we get the system of four 
equations. These equations can be solved considering the second terms in 
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the equations (23,24) as a small perturbation. In zero order we neglect by 
the second terms . The zero order solution has the form: 



jT £ Ue-^(s 2 +e 2 p) 



- \ qgiz(q,k z , u) - si#i q (q,£; 2 , u) 



J-oo A [ " 2 ' ^ kl-s\ 

. Jl9M,K,u) + s 2 g 2 >(c{,k z ,uj') . 
+s 2 (£ip-si) fc2 - g _ 2 (25) 



n f°° V j ipd / g^(q,A; 2 ,cu) - Sl ^ lq (q,A; 2 ,w) 



J.2 _ „2 

«, 2 6 X 



. ,q'92z{<l,k z ,u') + s 2 g 2 ,{q,k z ,uj') 

-s 2 (e lP + si) ri-^ f (26) 

fc 2 — s 2 



+cjs 2 ( Sl -p) — fc2 _ s _ 2 — >dk z , (27) 



P ^ J ipd/ - x#in(qA,^) 

-.'s 2 ( Sl + P ) 9 -^^ ^dk z , (28) 



«. = -^, «/. = ^, (29) 
p p 

where we have introduce the notation: 

A = e ^ d (e lP - Sl )(e 2 p - s 2 ) - e" tpd (e lP + Sl )(e 2 p + s 2 ), 

A' =e ipd ( Sl - p)(s 2 - p) - e- ipd (p + Sl )(p + s 2 ). 
The first order solution has the form: 

„■ = OflZfOA, K = _fel + f!)A, (30) 



where 



A 



A' = 



q y p 2 V 

uiq 2 
u>'q y V 
c 2 q 



(pe 2 - s 2 ) V y pd + (pej + sMe-** 



3 Calculation of the force of friction 

The frictional stresses a and —a which act on the surfaces of the two bodies 
can be obtained from the xz-component of the Maxwell stress tensor a^, 
evaluated at z = 0: 

i f + OO 

* = 8W-00 duj[{E3zE ^ + ^3^3.) + (B 3z B* 3x ) + (B* 3z B 3x )] z=0 (32) 

Here (..) denote a statistical average over the random field. The averaging 
is carried out with the aid of equation(4). Note that the components of the 
random field gi and g 2 referring to different media are statistically indepen- 
dent, so the average of their product gives zero. Writing the squares of the 
integrals (11-12) in the usual way as double integrals, and carrying out one 
integration over the 5— function, we obtain 

a = ^/^9K S 3z(q,w)^(q.w)) + <^q»w)S 3B (q,w)) 

+(B 3z (q,u;)B* 3x (q,u;) + (B* 3z (q, u)B 3x (q, u))] g ^>, (33) 

where one must substitute in place of E3 and B3 ,the expressions in the inte- 
grands (11-12) determined by the formulas (25-31), and the average product 
(gi(k,uj)gl(k,uj)} is to be taken as (1/2 + n(uj))e"(uj)5ik/TT 3 . For a given value 
of q it is convenient to express the components E x and B x in terms of the 
components along the two vectors e q and e n 

E x = (q x /q)E cl - (q y /q)E n , 
B x = (q x /q)B cl - (q y /q)B n . 

Thus we can write 
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+(5,(q,u;)E q (q,u;)) + (^(q, w)S q (q, w))] z=0 

[(^(q, W )K(q,w)> + <£;<?, w)£ n (q,")> 
+(S,(q,w)S;(q,c < ;)) + (S;(q, w)S n (q, w))] z=0 } . (34) 
Accordingly t) eq. (11-12) 



E z = 


(y z + Wz ) = (q/p)(w tl - u q ) = (qp*/ | p | 2 )(w q 


= 


^q + Wq, 


E n = 


+ Wn, 


B z = 


(cg/cj)(w n + u>„), 


= 


(cp/u)(w n - v n ), 


B n = 


(Ujp*/C |p| 2 )K - W q ). 



(35) 

After substituting these expressions into formula (34) one can see that 
the second term with q y is identically equal to zero . From equations (25-31) 
it follows that the zero and first order solutions are statistically independent, 
then, neglecting by the terms of the order (f/c) 2 ,from (34-35) we obtain 

o = lJ^ d ,J d \q x (\^P + P%\<\ 2 )-{\v^)) 



1\p\ 
+(P 



+ 



c\ 2 

CO J 



ip + p*)((\ w 



.0 |2 



- { V. 



|2 



)-(p-p*)<(«-<m ' 



(36) 

where we change the integration over dco between the limits — oo and +oo 
on the integration only over positive values of co what gives the extra factor 
two. 

Taking into account that p = p* for q < co/c , p = — p* for q > co/c, and 
carrying out the integration over dk z with the help of the formula 

dk 7 



I 

j — i 



,2|2 



'-oo \k\ — s 

after substituting (25-31) into (36) we obtain 

% 



s *y 



a 



8tt 3 



/ du d 2 qq x 

JO Jq<ui/c 
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(1- 




Rip 


2 )(l- 


R2p 


2 ) 




1 - e 2 ^ d R lp R 2p 


2 



x { i i .wo o- V - & - n( W )) + - i? s 

9 rt 2p I 



27T 3 JO Jq>uj/c 



lmRi p lmR 2p 



X [ I 1 - e-^'iM^, | 2 K ^ " &V) " ^ + [Rp ^ Rs] t (37) 



where 

_ £iP - Si _ Sj-p 

r^ip , j -Wis , i 

£j£> + Sj + p 

± _ efp - sf ± _ gf-p 
— ~i — ; — s"j — t - ; — 5 

£iP + Si Si +p 

ef{oo) = Ei{uo ± q x V) and sf (a;) = Si(aj ± q x V), % — 1, 2. Note that and 
.Ris are the electromagnetic reflection factors for p-polarized and s-polarized 
light, respectively, (p-polarized light has the electric field vector in the plane 
of incidence while the electric field vector is perpendicular to this plane for 
s-polarized light.) The first term in (37) is the contribution to the friction 
force from the propagating (radiating) electromagnetic field, i.e., the black 
body radiation. This term includes only the thermal radiation and is equal 
to zero at T = 0.. The second term is derived from the evanescent field, i.e., 
from the component of the electromagnetic field which decay exponentially 
with the distances away from the surfaces of the bodies. This term does not 
vanish even at T = OK because of quantum fluctuations in the charge density 
in the solids. 



4 Some limiting cases 

Let us first consider distances d « c/uj p , where uj p is the plasma frequency 
of the metals. For typical metals, c/uj p 200 A. In this case the main 
contribution comes from q » u p /c, and we have Si ~ s 2 ~ p ~ iq, R s ~ 
and ^ 

Rn ~ 

P E+l 

In this approximation the integration over d 2 q can be extended on the all q— 
plane. Using these approximations, the second term in (11) can be written 

as 

h 



47T 3 



duo J d 2 qq x e- 2Md 
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x { ( i i-T^hS5 2 - p + (1 " 2) ) - q * v) - n{uj)) + -> 



ft 

2^3 



dgj, J dq x q x e~ 2qd | y cL>[n(cj) - + g x i>)] 



x 



Imi?^Imi?2p 



e-2|p|di2+i22p 



2 




^^^^F^ 1 -')) « 38 » 

where we have used the relation n(—u) = —n(u) — 1. At zero temperature 
the Bose-Einstein factor n(cj) = and only the second term in (38) will con- 
tribute to the sliding friction; in this limit our expression for the friction force 
for two identical solids is in agreement with Pendry [his result is, however, 
not symmetric with respect to (1 <-> 2) which must be the case because of 
symmetry in the nonretarded limit]. In appendix A we show that the zero- 
temperature result can be generalized to include nonlocal optics effects, by 
replacing the reflection factor R p (ou) in (38) by the surface response function 
g(q,u). Next, let us consider the limiting case of low sliding velocity or high 
temperature namely, V « cd/dw, where dw = ch/ksT is the Wien length 
(typically d w ps 10 5 A). In this case 

T dn e ^/k B T hqx y 

n{u) - n{uj + q x V) « ~q x V— = — TT2 T^F 

and in the second term in (38) we can put 

n{uj) w k B T/huj 
Substituting these results in (38) gives 



hV f°° , , __ 2qd f°° ( _dn\ lmR lp lmR 2p 



POO POO / 

„■ / dqq 3 e- 2qd du - , . 
2W0 Jo \ dw) 1 - e~ 2 i d R^R 



lp-fl-2p 
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+ —k B T dq y / dg^e 9 / — — 39 

7r d j-oo io jo a; 1 — e lqa R\r,Ri n r 



2 



The second term in this expression is proportional to ~ V as V — > (see 
below) and can be neglected in the limit of small V. The first term is ~ V 
and is in agreement with the result obtained by Persson and Zhang if one 
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assumes local optics (which implies replacing g(q,u) — > R p {oj) in j2|). For 
free-electron like metals the local optics is accurate if d » I, where I is the 
electron mean free path in the metal. If this condition is not satisfied, the 
general formula of Persson and Zhang must be used. 

Let us consider two identical metals described by the dielectric function 

= 1 - , "}.. y (40) 

where r is the relaxation time and u p the plasma frequency. Thus, for small 
frequencies 

ImiL « — — ReR p « 1. 

UJ Z T 

p 

and if we neglect the imaginary part of R p in the denominator of the integrand 
in (13) we obtain 

cT = Jk B T\ 2 1 hV | 2n k B T 1 W 2 
\huOp J (uj p t) 2 d^ 45 huo p (uj p t) 2 u p d 5 ' 



where 



1 r°° dxx 2 

£=- ~ 0.5986. 

S 8io 



In deriving (41) we have used the following standard integrals 

r °° dxX 7T 2 r 00 dxX 3 7T 4 



6 jo c x — 1 15 

The ratio between the second and first term in (41) equals ~ (V/c)(dw/d), 
and in deriving (41) we have assumed that this quantity is much smaller than 
unity. As an example, if d = 10 A and V = lm/s then for typical metals at 
room temperature (ksT ~ 0.025 eV, uj p t w 100, huj p w 10 eV) the first and 
the second terms in (41) give a ~ 10~ 8 and w 10~ 13 N/m 2 , respectively. 
On the other hand, if (V/c)(dw/d) >> 1 we get 

UkBTV 1 nV +^ % Z(^\, (42) 



where 



huj p J (oo p t) 2 d A (u p t) 2 d A \doj P/ 
dxx 4 



2 9 7T 2 



0.024610. 
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The ratio of the second and first term in (42) equals ~ 0.l(V/c) 2 (d w /d) 2 . It 
is clear that at low temperature or high velocities, the second term in (42) 
will dominate. 

Next, let us consider the sliding friction to linear order in the sliding 
velocity when c/uj p « d « dw- There will be two contributions associated 
with R p and R s . As shown in Appendix B, the contribution from R p is 

3£hV ( d \ 2 k B T 1/1 d w \ 



7r 2 d 4 \dw J huj p uj p t \ e d 

The contribution a s from the term involving R s is given by (see Appendix 
B) 

<j s ps Cu p t—-^-, (44) 
c d z d w 

where C ps 0.394. Comparing (43) with (44) we obtain 

o s jop ps (u p T) 2 (hup/k B T) 2 . 

For typical metals at room temperature, Tiujp/kBT ~ 10 3 and uj p t ~ 100 so 
that <J s jo p ~ 10 10 , i.e., the main contribution comes from the term involving 
R s . As an illustration, if d = 10 4 A and V = lm/s then for metals at room 
temperature, characterized by the same parameter values as used above, one 
get a s PS 10~ 8 N/m 2 . 

Now, let us consider the radiative contribution to the friction force, which 
is given by the first term in (38). In linear order in the sliding velocity we 
get 



(Jrad = ^rC / du I d l qq 



f / " 1-1 X 

57T° JO Jq<w/c 



x i iiZ^w {-^l^^^l (45) 

For separations d much smaller than the Wien wavelength dw = ch/ksT 
we can put exp(ipd) ~ 1. In this case and for the small frequencies, when 
lu < k B T/h « l/r,we get for identical metals described by the dielectric 
function (40) 

1 (e*s) 2 + (es*) 2 1 e* + e 

I _^ r^j 

\ T~\ io ~ 77 ~r 



(1- 


Rp 






1 - 


-R 2 P 


2 



2 ' 4 | es | 2 2 I 

1 / UT\ 1 
-|H PS - 

2 V 2 7 2 



1 ^ ■ w n - i (46) 
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and the same result is obtained when R p in (46) is replaced by R s . The final 
result for the radiative friction force has the form . 

Grad = J^L [°° duouoM^) = —hv( 4 . (47) 
8n 2 c 4 Jo V ; 120 V hc J 

Note that the radiative stress does not depend on the separation and is pro- 
portional to T 4 .The latter result is, of course, only valid as long as d is small 
compared with the lateral extend (or linear size) L of the bodies. When d 
becomes comparable with or larger than L, the friction force between the two 
bodies will decrease monotonically with increasing d. At room temperature 
and at the sliding velocity V = lm/s one get a rad pa 10 _15 N/m 2 . The ratio 
of this contribution to a s from (44) is 

<J S \d\V J flUJp LdpT \dw 

Thus for d ~ dw (which is of order ~ 10 5 for typical metals at room tem- 
perature) the nonradiative part dominate over the radiative contribution by 
a factor ~ 10 6 . However, for large enough distances the radiative part dom- 
inate as this contribution is finite for arbitrary separations d. 




5 Summary and conclusion 

We have calculated the friction force between two arbitrary bodies with flat 
surfaces separated by a vacuum slab of thickness d, and moving with a rel- 
ative velocity V. The separation d is assumed to be so large that the only 
interaction between the bodies is via the electromagnetic field associated 
with thermal or quantum fluctuations in the solids. A general formula for 
the friction force has been obtained, which is valid for arbitrary velocity 
V, separation d and temperature T, and applicable to any bodies. At low 
sliding velocity only thermal fluctuations give a contribution to the friction 
force, linearly proportional to the velocity V. Quantum fluctuations give a 
nonlinear (in V) contribution to the friction force and is usually negligible 
compared with the thermal contribution. [There is also a contribution from 
quantum fluctuations which is proportional to V, resulting from higher order 
electron-photon processes than considered in our work. However, this con- 
tribution decays as ~ exp(—2Gd) (where G = 2n/a is the smallest reciprocal 
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lattice vector) and is negligible small already for d = 10 A. See [TT| |2] .] We 
have studied the detailed distance dependence of the friction force from short 
distances, where retardation effects can be neglected, to large distances where 
retardation effects and black body radiation are important. In most practical 
cases, involving sliding of a block on a substrate, the van der Waals friction 
makes a negligible contribution to the friction force (the main part of the fric- 
tion arises from the regions of real contact between the solids). However, in 
some special cases the van der Waals friction is very important |2|. For exam- 
ple, quantum fluctuations contribute in an important manner to the friction 
force acting on thin physisorbed layers of atoms sliding on metallic surfaces 
[T2] [In this case there is an overlap of the wavefunctions of the sliding layer 
and those of the metal, which result in second contribution to the friction 
force, derived from the repulsive "contact interaction" (Pauli repulsion) be- 
tween the sliding layer and the substrate] . In addition, the contribution from 
thermal fluctuations gives the dominating drag force in some experiments 
involving parallel 2D-electron systems IT51 
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In this appendix we derive an expression for the (nonlinear) sliding friction 
using a nonlocal optics description of the metals. For simplicity, we focus 
on zero temperature T = OK and assume that d is so short that retardation 
effects can be neglected (only in this "short-distance" region will nonlocal 
effects be important). The calculation is based on the formalism developed in 
Let us first define the linear response function g(q, u) which is needed 
below. Assume that a semi- infinite metal occupy the half space z < 0. A 
charge distribution in the half space z > d gives rise to an (external) potential 
which must satisfy Laplace equation for z < d and which therefore can be 
written as a sum of evanescent plane waves of the form 

<P ext = foe^e^*-^ 

where q = {q x ,Q y ) is a 2D-wavevector. This potential will induce a charge 
distribution in the solid (occupying z < 0) which in turn gives rise to an 
electric potential which must satisfy the Laplace equation for z > 0, and 
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which therefore can be expanded into evanescent plane waves which decay 
with increasing z > 0. Thus the total potential for < z < d can be 
expanded in functions of the form 



where the reflection factor g = g(q,u). For the present purposes, we can 
treat the low-energy electron-hole pair excitations in the metals as bosons. 
As shown in Ref. the Hamiltonian for the total system can be written 



as 



H = E ^Vikq Ql frqai + E ^ 'qo^qaa^Va + ^Ujb + b 
qctl qc*2 

+ E C qai e-^ (& qQl e^ x " +vt ) + h.c.) . (48) 

Here u; qQl , 6+ and fe qai are the angular frequency and creation and annihila- 
tion operators for the bosons (of solid 1) with the quantum numbers (q,ofi), 
and C qQl parameters determining the coupling between the boson excitations 
in solid 1 with the electrons in solid 2. Similarly, b+ a2 and 6 qQ2 are creation 
and annihilation operators for bosons in solid 2, and (x n , z n ) is the position 
operator of electron n in solid 2, which in principle could be expressed in 
terms of the operators and 5 qQ , 2 , but for the present purpose this is not 
necessary. As shown in [14j, C qCfl can be related to lmgi(q, uS) via 

2e 2 h 

E I C n<*i I 2 5 ( u ~ = —7 — Im^i(g, u) (49) 
We can write the interaction Hamiltonian between solid 1 and 2 as 

H' = E ( V 1 ™ + h.c) 
q 

Using time-dependent perturbation theory (with H as the perturbation) 
we can calculate the energy transfer from the translational motion (kinetic 
energy) to internal excitations in the solids (boson excitations u; qai and u; qQ2 
in solid 1 and 2, respectively): 

2ir 2 

" qai«2 

xe' ¥ | (n qQl = 1, n qa2 = 1 | E e'^^e-^b^ | 0, 0) | 2 (50) 
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where u; q =| q • V |. To simplify (A3), let us write 

5(w q + ^ qQ2 - u mi ) = J du'5(uj' - w qQl )<5(u; q - to' - tu qQ2 ) (51) 
Substituting (A4) in (A3) and using (A2) gives 

p = W ^ u^ e _ 2qd f du , lmgi{q ^ ^ )MgK _ ^ (52) 
A q g J 

where 

M» = £ % - ^ 2 ) I = 1 I E e^^e^ | 0) | 2 

Q2 It 

But it has been shown elsewhere that [HJ 

-l^Im^(g^) = I K* = 1 I E e ^""" e " iq,X " I °) I' 



M,H = — ^Im^(g,^) (53) 



so that 

2?r 2 e 

Substituting this result in (A5) gives 

P = — E ^ e ' 2Qd / ^' Im 3i(g, u')Img 2 (q, ^ - u') (54) 

7T q J 

Finally replacing 

?^/*» 

and using the relation P = aAV between the power P and the shear stress 
a gives 



2tt 3 



J d 2 q | q x | e" 29d y o ' cL/Imc^, u )1mg 2 (q, \ q x \ V - J) (55) 



where we have used that lmg(q, u) = for uo < O.The coupling H' does not 
only give rise to real excitations but also to screening (image charge effects). 
To take these into account one must go to higher order in perturbation theory. 
Following Ref.j2] this gives a modification of (A8) 

h f 2 | q x | e- 2qd 

a = —r a q 



2tt 3 J *| 1- g 1 (q 1 0)g 2 (q,0)e- 2 i d | 2 
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r\qx\v 

x / du'lmg 1 (q,uj')lmg2(q, \ q x \ V - u) (56) 
Jo 

where we have assumed that the small frequencies involved in the real excita- 
tions are screened in an adiabatic manner so that g\ and (72 can be evaluated 
at zero frequency in the screening factor. In the case of local optics, this 
expression for a agree with the last term in (38) evaluated at zero tempera- 
ture. At finite temperature (T > OK) an extra factor of [1 + 2n{uj )] must be 
inserted in the frequency integral in (A9) to take into account the enhanced 
probability for excitation of bosons at finite temperature. For T > OK one 
must, in addition to the process considered above, also include scattering 
processes where a thermally excited boson is annihilate either in solid 1 or 
in solid 2, namely (n qcri = 0,n qQ2 = 1) — > (1,0) and (1,0) — > (0, 1). These 
processes was considered in Ref. and give, in the local optics case, the 
frictional stress corresponding to the first term in (38). 

In this appendix we calculate the sliding friction to linear order in the 
sliding velocity when c/uj p « d « d w (where dw = ch/ksT). In this case 
the main contribution a comes from the first term in (39) to which we must 
add the similar term involving R s . In this integral we replace the integration 
variable q with p = 2dq. The integral over p is divided into two parts: the 
integral over (0,po) and over (p Q ,oo), where po ~ dkBT/ch « 1. In the 
integral over (p , 00) and for u > u , where 



2 1 



^0 

we can put 



u) v d) r 




(58) 



uj p d \ujtJ 



Let us first consider the contribution a v to (39) from terms involving R p . We 
then obtain 

a P = n -A-^ r^{-^)r^{— V 

p 7T 2 \2d) Jpo FF dp\ eP-lJJo 2uo 2 p t\ c J 

3 hv / r 1 ,_ 1 r°° A (k B Td\ 2 k B T 1 



1 + / dp^— + / dpe~ p 11 1 , 

87r z ct 4 V Jp e p — 1 Ji / \ nc J nujp uj v t Jo 
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SZhv fk B Td\ 2 k B T 1 / | 1 | ^d w \ 

7T 2 rf 4 \ flC J hujp UJpT \ e d J 

where £ = 0.5986 (see Sec. 3). The integral over < p < p can be shown 
to give a negligible contribution to the linear (in the sliding velocity) friction 
force. This follows from the following equation 

d (lmR p ) 2 _ d (lmR p ) 2 _ d ( (lm(s/e)f\ Q 



duj | 1 — e p R p Rp | 2 du \ 1 — -R p i? p | 2 da; \ | s/e | 2 / 

where we have used that Im(s/e)/ | s/e \ is approximately independent of 
frequency. Next, let us consider the contribution a s to (39) from the term 
involving R s . We get 



2 1 

5 Tiv % ( c \ f°° dpp A f°° dx d ( 1 



167T 2 d 4 ksTr \cu p d J Jo e p — 1 J xo x du 
hv ( ksTT\ I uo p d\ 2 



C 



d 4 \ h 



where 



5 

32tt 2 Jo eP-1 

hujQ h ( c \ 2 

k B T k B Tr \uj p dj 
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